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The theory of scale relativity provides a new insight into the origin of fundamental laws in 
physics. Its application to microphysics allows to recover quantum mechanics as mechanics on a 
non-differentiable (fractal) space-time. The Schrodinger and Klein-Gordon equations have already 
been demonstrated as geodesic equations in this framework. We propose here a new development 
of the intrinsic properties of this theory to obtain, using the mathematical tool of Hamilton’s bi¬ 
quaternions, a derivation of the Dirac equation, which, in standard physics, is merely postulated. 
The bi-quaternionic nature of the Dirac spinor is obtained by adding to the differential (proper) 
time symmetry breaking, which yields the complex form of the wave-function in the Schrodinger 
and Klein-Gordon equations, the breaking of further symmetries, namely, the differential coordinate 
symmetry (dx^ e-> —dx^) and the parity and time reversal symmetries. 

PACS numbers: 03.65.Pm, 12.90.-|-b 


I. INTRODUCTION 

The theory of scale relativity consists of generalizing to scale transformations the principle of relativity, which has 
been applied by Einstein to motion laws. It is based on the giving up of the assumption of space-time coordinate 
differentiability, which is usually retained as an implicit hypothesis in current physics. Even though this hypothesis 
can be considered as roughly valid in the classical domain, it is clearly broken by the quantum mechanical behavior. 
It has indeed been pointed out by Feynman (see, e.g., Ref. |Q|) that the typical paths of quantum mechanics are 
continuous but non-differentiable. 

In the present paper, after a reminder about the Schrodinger and Klein-Gordon equations, we apply the scale- 
relativistic approach to the issue of the understanding of the nature of bi-spinors and of the Dirac equation in a 
space-time representation. A step in this direction has been made by the work of Gaveau et al. Q who have general¬ 
ized Nelson’s stochastic mechanics to the relativistic case in (1-1-1) dimensions. However, an analytic continuation 
was still needed to obtain the Dirac equation in such a framework, and, furthermore, stochastic mechanics is now 
known to be in contradiction with standard quantum mechanics 0. 

Basing himself on the concept of fractal space-time as a geometric analog of quantum mechanics, Ord has de¬ 
veloped a generalized version of the Feynman chessboard model which allows to recover the Dirac equation in (3-1-1) 
dimensions, without analytical continuation @-0- We develop in the present paper an approach involving also a 
fractal space-time, but in a different way: namely, the fractality of space-time is derived from the giving up of its 
differentiability, it is constrained by the principle of scale relativity, and the Dirac equation is derived as an integral 
of the geodesic equation. This is not a stochastic approach (the use of statistically defined quantities is derived and 
we do not use Fokker-Planck equations), so that it does not come under the contradictions encountered by stochastic 
mechanics. 

One of the fundamental reasons for jumping to a scale-relativistic description of nature is to increase the general¬ 
ity of the first principles retained to construct the theory. Since more than three centuries, i.e., since Newton and 
Leibnitz, the founders of the integro-differential calculus, physics relies on the assumption that coordinates, and more 
generally all the physical equations, are a priori differentiable. Indeed, one of the most powerful avenue for reaching a 
genuine understanding of the laws of nature has always been to construct them by attempting to lessen the number of 
a priori assumptions retained, a principle known as Occam’s razor. From this point of view, the laws and structures 
of nature are simply the most general laws and structures that are physically possible. It culminates nowadays in 
Einstein’s principle of “general” relativity and his description of gravitation from its analogy with the manifestation 
of a Riemannian geometry of space-time, obtained by giving up the Euclidean flatness hypothesis. 

However, in the light of the above remark, Einstein’s principle of relativity is not yet fully general, since it applies to 
coordinate transformations that are continuous and, at least two times, differentiable. The aim of the theory of scale 
relativity is to bring to light laws and structures that would be the manifestation of still more general transformations, 
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namely, continuous ones, either differentiable or not. The standard “general relativistic” theory will be automatically 
recovered as the special differentiable case, while new effects would be expected from the non-differentiable part. 

Giving up the assumption of differentiability has important physical consequences: one can show Bi that spaces 
of topological dimension Dt, which are continuous but non-differentiable, are characterized by a Dt measure which 
becomes explicitly dependent on the resolution (i.e., the observation scale) e at which it is considered and tends to 
infinity when the resolution interval e tends to zero. Therefore, a non-differentiable space-time continuum is necessar¬ 
ily fractal, in the general meaning given to this word by Mandelbrot |p^ , who thus means, not only self-similar, but 
explicitly scale-dependent in the more general way. This statement naturally leads to the proposal of a geometric tool 
adapted to construct a theory based on such premises, namely, a fractal space-time (see, e.g.. Refs. [^, [ll| , |l^ -p^ ). 
This allows to include the resolutions in the definition of the state of the reference system and to require scale co- 
variance of the equations of physics under scale transformations. A main consequence is that the geodesics of such 
a non-differentiable space-time are themselves fractal and in infinite number between any two points. This strongly 
suggests that the description could borrow some tools from the mathematical formalism of statistics, i.e., will provide 
a probabilistic interpretation, therefore naturally extending the realm of quantum behaviour to a larger spectrum of 
available scales. 


The domains of application are typically the asymptotic scale ranges of physics: (i) Microphysics: small length-scales 
and time-scales, (ii) Cosmology: large length-scales. In this framework, new solutions can be brought to the problems 
of the vacuum energy density, of Mach’s principle and of Dirac’s large number hypothesis, leading to a theoretical 
prediction of the value of the cosmological constant |17|. (iii) Complex chaotic systems: long time-scales. It has been 
suggested that a general theory of structure formation could take a quantum-like form, but with an interpretation 
differing from that of standard quantum mechanics 00- Several applications to the problem of the formation and 


extra-solar plan¬ 
binary galaxies M 


evolution of gravitational structures have been completed in the recent years: solar system 111 

etary systems planets around ^Isars |Q , satellites of giant planets |23| , binary stars ^ ^ ^, 

large-scale distribution of galaxies 1^]. Moreover, an additional scalar field appears naturally in such a construction 
(as a manifestation of the fractality of space-time, in analogy with the Newton field arising from its curvature), which 
could explain the various effects usually attributed to unseen, “dark matter” . It has also been suggested that this 
approach could be applied to the morphogenesis of a wide class of hierarchical structures, even outside the strict field 
of physics: for instance, in biology, structures are “density waves” rather than solid objects and could therefore be 
described in terms of probability amplitudes, being themselves solutions of a generalized Schrodinger equation . 


The theory of scale relativity can also be considered, from the point of view of its construction, as a completion of 
the standard laws of classical physics (motion in space, displacement in space-time) by new scale laws, in which the 
space-time resolutions are used as intrinsic variables, playing for scale transformations the same role as velocities for 
motion transformations. We hope such a stage of the theory to be only provisional, and the motion and scale laws to 
be treated on the same footing in the complete theory. But before reaching such a goal, various possible combinations 
of scale and motion laws, leading to a number of sub-sets of the theory, have to be developed. 


We can first distinguish: (1) Pure scale laws: description of the internal structures of a non-differentiable space-time 
at a given point. (2) Induced effects of scale laws on the equations of motion: generation of the quantum mechanichal 
behavior as motion in a non-differentiable space/space-time. (3) Scale-motion coupling: effects of dilations induced 
by displacements, that have been tentatively interpreted as gauge fields, allowing to suggest the existence of new 
relations between the masses and the charges of elementary particles 0: 


Several levels of the description of pure scale laws (1) can also be distinguished, which therefore induce new sub-sets 
inside each of the three above cases. These levels can actually be considered as reproducing the historical development 
of the theory of motion: (a) Galilean scale relativity: laws of dilation exhibiting the structure of a Galileo group, i.e., a 
fractal power law with a constant fractal dimension, (b) Special scale relativity: generalization of the laws of dilation 
to a Lorentzian form. The fractal dimension becomes a variable and is assimilated to a fifth dimension. A finite 
impassable length/time-scale, invariant under dilations, shows up. It owns all the physical properties of a zero (at 
small scales) or of an infinite (at large scales) - an infinite energy-momentum would be needed to reach it - and it plays 
for scale laws the same role as played by light velocity for motion, (c) Scale dynamics: while the first two above cases 
correspond to “scale freedom”, one can also consider the distorsion from strict self-similarity which would proceed from 
the effect of a “scale force” or “scale field” [^,^. (d) General scale relativity: in analogy with the gravitation field 
being assimilated to the effect of the space-time geometry in Einstein’s general motion relativity, a future description 
of the scale field of (c) could be done in terms of a Riemannian geometry of the five-dimensional scale space, therefore 
assuming differentiability of scale transformations, (e) Quantum scale relativity: the scale transformations are still 
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assumed to be continuous but no more differentiable (as it has been previously done for space-time). One is therefore 
confronted, for scale laws, to the same conditions that lead to a quantum-like behaviour in a fractal space-time, and 
may conjecture that quantum scale laws could emerge at this (to be developed) level of the theory. 

A forth spliting of the different levels and sub-levels described above is to be considered, depending on the status 
of the motion laws: non-relativistic, special-relativistic, general-relativistic. 


In the present work, we focus our attention on the microphysical scale of motion in a non-differentiable space-time, 
in the framework of Galilean scale relativity (i2a). Therefore, we recover, from the first principles of this theory, the 
three main evolution equations of standard quantum physics as geodesic equations on a fractal space/space-time: 

- the Schrodinger equation 

- the (free) Klein-Gordon equation 

- the (free) Dirac equation 

which are here demonstrated, while they are merely postulated in standard quantum mechanics (the compatibility 
between the Dirac equation and the scale-relativistic formalism has already been studied by Pissondes j^). 


A reminder of the derivation of the Schrodinger equation is given in Sec. ||, where we actualize the former works 
dealing with this issue, proposing a more accurate and profound interpretation of the nature of the transition from 
the non-differentiable (small scales) to the differentiable (large scales) domain, and carefully justifying the different 
key choices made at each main step of the reasoning. An analogous reminder and updating is pro posed in Sec. |I|, for 
the free Klein-Gordon equation exhibiting a complex (number) form. Then we derive, in Sec. IV, the bi-quaternionic 
form of the same Klein-Gordon equation, from which the Dirac equation naturally proceeds, as shown in Sec. The 
two appendices are devoted to the presentation and justiffcation of the mathematical tools used in this framework. 


The reader already familiar with the scale-relativistic formalism, and mainly int erest ed in the broad new results 
presented here, can limit his attention to Secs. 0 and 0, after a detour via Sec. |[I C| , where a new symbolism is 
introduced. 


II. THE SCHRODINGER EQUATION REVISITED 

In previous works, the Schrodinger equation has already been derived, in the framework of scale relativity, as a 
geodesic equation in a fractal three-space (see, e.g.. Refs. ]ll] , p^ , p^ ). We give, in the present section, a reminder of 
the main steps of the reasoning leading to its formulation, therefore making the reader familiar with the structure 
and tools of the theory. We also update and give a more precise meaning to some of the most subtle issues we are 
dealing with. 


A. Scale invariance and Galilean scale relativity 


A power law scale dependence, frequently encountered in natural systems, is described geometrically in terms of 
fractals ]r^ , ^l[ |, and algebrically in terms of the renormalization group 
invariant laws can be identiffed with the “Galilean” version of scale-relativistic laws. 


331. As we show below, such simple scale- 


Consider a non-differentiable (fractal) curvilinear coordinate C(x^e), that depends on some space-time variables x 
and on the resolution e. Such a coordinate generalizes to non-differentiable and fractal space-times the concept of 
curvilinear coordinates introduced for curved Riemannian space-times in Einstein’s general relativity [ p]| . Rather 
than considering only the strict non-differentiable mathematical object C{x), we are interested in its various reso¬ 
lutions e. Such a point of view is particularly well-adapted to applications in physics, as any real measurement is 
always performed at a finite resolution. In this framework, C{x) becomes the limit when e —> 0 of the family of 
functions C{x, e). But while C{x, 0) is non-differentiable, C{x, e), which we call a fractal function, is differentiable for 
all e yf 0 (see Refs. for a precise mathematical definition of a fractal function). The physics of a given process 

will therefore be completely described if we succeed in knowing L{x, e) for all values of e. jC(x, e), being differentiable 
when e yf 0, can be the solution to differential equations involving the derivatives of £ with respect to both x and e. 

Let us now apply an infinitesimal dilation e ^ = e(l -I- dp) to the resolution. Being, at this stage, interested in 

pure scale laws, we omit the x dependence in order to simplify the notation and obtain, at first order. 


3 






C{t') = C{e + edp) = £(e) + 


dCje) 

de 


edp = (1 + Ddp)C(e), 


( 2 . 1 ) 


where D is, by definition, the dilation operator. The identification of the two last members of this equation yields 

( 2 . 2 ) 


D= ^ ^ 

^ de 9 In e ' 


This well-known form of the infinitesimal dilation operator shows that the “natural” variable for the resolution is 
Ine, and that the expected new differential equations will indeed involve quantities as dL{x,e)/dhie. The equations 
describing scale dependence are known to be the renormalization group equations, as first proposed by Wilson 


We intend to limit the present approach to the consideration of the simplest form that can be exhibited by this 
class of equations, leaving to future works the study of more complete cases. We therefore interpret the simplest 
renormalization group-like equation implying an essential physical quantity like C as stating that the variation of C 
under an infinitesimal scale transformation d In e depends only on C itself; namely, C determines the whole physical 
behavior, including the behavior under scale transformations. We thus write 


dC{x, e) 
9 In e 


= m- 


(2.3) 


Still interested in the simplest form for such an equation, we expand (3{C) in powers of £. This can always be done 
since C may be renormalized, dividing it by its largest value, in such a way that the new variable C remains <C 1 in 
its variation domain. We obtain, to the first order, the linear equation 


of which the solution is 


dC{x, e) 
91ne 


C{x,e) = Cq{x) 


= a + bC 




-b 


(2.4) 


(2.5) 


where A °C,{x) is an integration constant and £o = —a/b. 


For b < 0, Eq. ( |2.5| ) gives a fractal scale-invariant behavior at small scales (e <C A), with a fractal dimension 
D = Dt — b. For a curvilinear coordinate £, the topological dimension Dt is equal to 1. The scale dimension, 
6 = D — Dt, is defined, following Mandelbrot as 


S = 


dlnC 


dln(A/e) 

When 6 is constant, one obtains an asymptotic power law resolution dependence 
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C{x,e) = Co{x) 


( 2 . 6 ) 


(2.7) 


We have stressed in the introduction the fact that the new scale laws to be constructed were expected to come 
under the principle of relativity extended to scale transformations of the resolutions e. Let us now check that this is 
indeed the case for such a simple self-similar scaling law. The Galilean structure of the group of scale transformations 
which correspond to this law can be verified, in a straightforward manner, from the fact that the involved quantities 
transform, under a scale transformation e —> e', as 


C{e') C{e) e 

In = In ^ + (5(e) In - , 
£o £o e 


( 2 . 8 ) 


5{e')=5{e). 


(2.9) 


These transformations have exactly the structure of the Galileo group, as confirmed by the dilation composition 


law. 


e”, which writes 
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e e 


1 


(2.10) 


It is worth noting that Eq. (2.5) gives also a transition from a fractal to a non-fractal behavior at scales larger than 
some transition scale A. 


The solutions corresponding to the case 5 > 0 yield scale dependence at large scales, broken to scale independence 
below A. While 5 < 0 is characteristic of the microphysical situation (it yields both quantum phenomena and the 
quantum-classical transition), the case 6 > 0 is encountered in the astrophysical and cosmological domains (see, e.g., 
Ref. @). 

We want to stress that we have re tain ed here the simplest case of a first order limited expansion of /3(>C) in powers 
of jC. Going to higher order in Eq. (2.4) would produce an increasing number of transition scales: two for a second 
order limited expansion, three for a third order, and so on. 


B. Lagrangian approach to scale laws 


We are naturally led, in the scale-relativistic approach, to reverse the definition and meaning of the variables. The 
scale dimension S becomes a fundamental variable. It plays, for scale laws, the same role as played by (proper) time in 
motion laws. It remains constant in only very peculiar situations, namely in the case of scale invariance, corresponding 
to “scale freedom”, which is the (a) case here studied. When going to a more general description of pure scale laws 
(b-c-d-e), one is led to consider a scale dimension 5(e), which is no more a constant [p|. 


The resolution, e, can therefore be defined as a derived quantity in terms of the fractal coordinate £ and of the 
scale dimension 5 


E = ln(^) 


dln£ 

dS 


( 2 . 11 ) 


Our identification of the standard fractal behavior with Galilean scale laws can now be fully justified. We assume 
that, as in the case of motion laws, the scale laws can be constructed from a Lagrangian approach. A scale Lagrange 
function £(ln£, V, S) is introduced, from which a scale action is constructed 


S = 


r^2 


L(ln £, V, 6)dS. 


( 2 . 12 ) 


The application of the action principle yields a scale Euler-Lagrange equation that writes 


d dL dL 
d^W ^ 51n£ ■ 


(2.13) 


The simplest possible form for the Lagrange function is a quadratic dependence on the “scale velocity”, V, (i.e., 
L oc tA^) and the absence of any “scale force” (i.e., dL/d\nC — 0), which is the equivalent for scale of what inertia is 
for motion. Note that this form of the Lagrange function becomes fully justified, as in the case of motion laws, when 
one jumps to special scale relativity, then back to the Galilean limit The Euler-Lagrange equation becomes in 
this case 


U 


= 0 y = const. 


(2.14) 


The constancy oi V = ln(A/e) means that it is independent of the “scale time” 6. Equation (2.11) can therefore 
be integrated to give the usual power law behavior, £ = £o(A/e)'^. This reversed viewpoint has several advantages 
which allow a full implementation of the principle of scale relativity: 

(i) The scale dimension 5 is giv en its actual status of “scale time” and the logarithm of the resolution, V, its status of 
“scale velocity” (see Eq. (2.11)). This is in accordance with its scale-relativistic definition, in which it characterizes 
the “state of scale” of the reference system, in the same way as the velocity v = dx/dt characterizes its state of motion. 

(ii) This leaves open the possibility of generalizing our formalism to the case of four independent space-time resolutions, 
= ln(A'"/e'") = dlnC>^/dS. 

(iii) Scale laws more general than the simplest self-similar ones can be derived from more general scale Lagrangians 
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C. Transition from non-differentiability (small scales) to differentiability (large scales) 


The aim of the present subsection is to identify the consequences of the giving up of the coordinate differentiability. 
We develop here the formalism relevant for the derivation of the non-relativistic motion Schrodinger equation, valid in 
a fractal three-space with the time t as a curvilinear parameter. This formalism will be easily generalized, in Sec. 01, 
to relativistic motion involving a four-dimensional space-time with the proper time s as a curvilinear parameter. 


Strictly, the non-differentiability of the coordinates implies that the velocity 

dX X(t + dt)-X(t) 

V = —— = lim -;- 

dt dt^o at 


(2.15) 


is undefined. This means that, when dt tends to zero, either the ratio dX/dt tends to infinity, or it fluctuates without 
reaching any limit. However, as recalled in the introduction, continuity and non-differentiability imply an explicit 
scale dependence of the various physical quantities, and therefore of the velocity, V. We thus replace the differential, 
dt, by a scale variable, 6t, and consider V(t, St) as an explicit fractal function of this variable. In the case of a constant 
fractal dimension D, the resolution in X, e{6t), corresponds to the resolution in t, St, according to the relation [pd]| 




(2.16) 


The advantage of this method is that, for any given value of the resolution, St, differentiability in t is recovered, 
which allows to use the differential calculus, even when dealing with non-differentiability. However, one should be cau¬ 
tious about the fact that the physical and the mathematical descriptions are not always coincident. Indeed, once St is 
given, one can write mathematical differential equations involving djdt, make dt —f 0, then solve for it and determine 
V(t, St). Actually, this is a purely mathematical description with no physical counterpart, since, for the system under 
study, the very consideration of an interval dt < St, as occuring in an actual measurement, changes the function V 
(such a behavior, described by Heisenberg’s uncertainty relations, is experimentally verified for any quantum system). 
However, there is a peculiar subspace of description where the physics and the mathematics coincide, namely, when 
making the particular choice dt = St. We work, for the Schrodinger and complex Klein-Gordon equations, with such 
an identihcation of the time differential and of the new time resolution variable, but are led to give it up in Secs. IV 
and 0. 

The scale dependence of the velocity suggests that we complete the standard equations of physics by new differential 
equations of scale. We therefore apply to the velocity a nd t o the differential element, now interpreted as a resolution, 
the reasoning applied to the fractal function C in Sec. HA. Writing the simplest possible equation for the variation 
of the velocity V{t, dt) in terms of the new scale variable dt, as a first order renormalization group-like differential 
equation of the form of Eq. (^.3|); then, Taylor-expanding it, as in Eq. (^^), using the fact that V < 1 (in motion- 
relativistic units c = 1), we obtain the solution 


V = V + w = V 




l-ljD 


(2.17) 


where we have set b = 1/D — 1. Here, v is the large-scale part of the velocity, LS < V >, (see below the definition of 
the large-scale part operator LS), w is an explicitly scale-dependent fractal fluctuation and r and C are chosen such 
that LS < C >= 0 and LS < >= 1. 


We recognize here the combination of a typical fractal behavior, with a fractal dimension D, and of a breaking 
of the scale symmetry at the scale transition t. As we shall see, in what follows, r will be identified with the de 
Broglie scale of the system (r = fi/E), since V ~ v, when dt t (classical behavior), and V ~ w, when dt ^ t 
(fractal behavior). Recalling that D = 2 plays the role of a critical dimension j^, we stress that, in the asymptotic 
scaling domain, w oc {dtfT)~^^^, in agreement with Ref. Q for quantum paths, which allows to identify the fractal 
domain with the quantum one. In the present paper, only these simplest scale laws with a fractal dimension D — 2 
are considered. As recalled in Sec. |H A| , these laws can be identihed with “Galilean” approximations of more general 
scale-relativistic laws where the fractal dimension becomes itself variable with scale. Therefore, all the results that 
will be obtained here (Schrodinger, Klein-Gordon and Dirac equations) are open to new generalizations, which we 
shall detail in forthcoming works. 

The above description strictly applies to an individual fractal trajectory. Now, one of the geometric consequences of 
the non-differentiability and of the subsequent fractal character of space/space-time itself (not only of the trajectories) 
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is that there is an infinity of fractal geodesics relating any couple of points of this fractal space. It has therefore been 
suggested [Q that the description of a quantum mechanical particle, including its property of wave-particle duality, 
could be reduced to the geometric properties of the set of fractal geodesics that corresponds to a given state of this 
“particle”. In such an interpretation, we do not have to endow the “particle” with internal properties such as mass, 
spin or charge, since the “particle” is not identified with a point mass which would follow the geodesics, but its internal 
properties can simply be defined as geometric properties of the fractal geodesics themselves. As a consequence, any 
measurement is interpreted as a sorting out of the geodesics of which the properties correspond to the resolution scale 
of the measuring device (as an example, if the “particle” has been observed at a given position with a given resolution, 
this means that the geodesics which pass through this domain have been selected) . 


The transition scale appearing in Eq. ( 2.17 
resolution at which it is considered (see Figs 
dX, of the system as a sum of two terms 


2.17) vie 


ields two dis tinct behaviors of the system (particle) depending on the 
Equation (2.17) multiplied by dt gives the elementary displacement. 


dX = dx + d^, 


(2.18) 


d^ representing the fractal fluctuations around the “classical” or “large-scale” value, dx, each term being defined as 


dx = V dt, 


(2.19) 


d^ = , ( 2 . 20 ) 

which becomes, for D = 2, 

dC = ( 2 . 21 ) 


with 227 = To = Tv'^. Owing to Eq. (2.17), we identify r as being the Einstein transition scale, Ti/E = and 

therefore, as we shall see further on, 227 = tq is a scalar quantity which can be identified with the Compton scale, 
fi/mc, i.e., it gives the mass of the particle up to fundamental constants. 


We note, from Eqs. ( 2.18 ) to ( 2.21 ), that dx scales as dt, while d^ scales as dt^. Therefore, the behavior of the 
system is dominated by the d^ term in the non-differentiable “small-scale” domain (below the transition scale), and 
by the dx one in the differentiable “large-scale” domain (above the transition scale). The logarithmic dependence of 
dx and d^ with respect to dt is illustrated in Fig. |^. 
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FIG. 2. Small-scale to large-scale transition, velocity 



FIG. 3. Small-scale to large-scale transition, fractal dimension 


Now, the Schrodinger, Klein-Gordon and Dirac equations give results applying to measurements performed on 
quantum objects, but realised with classical devices, in the differentiable “large-scale” domain. The microphysical 
scale at which the physical systems under study are considered induces the sorting out of a bundle of geodesics, 
corresponding to the scale of the systems (see above), while the measurement process implies a smoothing out of the 
geodesic bundle coupled to a transition from the non-differentiable “small-scale” to the differentiable “large-scale” 
domain. We are therefore led to define an operator LS < >, which we apply to the fractal variables or functions 







each time we are drawn to the “large-scale” domain where the dt behavior dominates. The effect of LS is to extract, 
from the fractal variables or functions to which it is applied the “large-scale part”, i.e., the pa rt sca ling as dt. This 
justify our above writting LS < ( >= 0 and LS < >= 1, as it is straightforward, from Eq. ( p.l7 ), that the terms 

involving ( scale asdt^/^ and those involving scale as dt. This also leads us to state, for the a dimensionless 
coefficient in Eq. (2.21), LS < a >= 0 and LS < >= 1. 


Please note the improvement of our new definition in terms of the large-scale part with respect to the previous 
interpretation in terms of an averaging process ||ll|] . 


D. Differential-time symmetry breaking 


Another consequence of the non-differentiable nature of space (space-time) is the breaking of local differential 
(proper) time reflection invariance. The derivative with respect to the time t of a differentiable function / can be 
written twofold 


# ^ f{t + dt) - f{t) ^ fit) - f{t - dt) 
dt dt^o dt dt^o dt 


( 2 . 22 ) 


The two definitions are equivalent in the differentiable case. One passes from one to the other by the transformation 
dt —dt (local differential time reflection invariance), which is therefore an implicit discrete symmetry of differentiable 
physics. In the non-differentiable situation, both definitions fail, since the limits are no longer defined. In the new 
framework of scale relativity, the physics is related to the behavior of the function during the “zoom” operation on 
the time resolution St, identified with the differential element dt. Two functions /(_ and f'_ are therefore defined as 
explicit functions of t and dt 


r'u f{t + dt,dt)-f{t,dt) 

f+[t,dt)- 

(2.23) 

1 

1 

II 

(2.24) 


When applied to the space coordinates, these definitions yield, in the non-differentiable domain, two veloc¬ 
ities that are fractal functions of the resolution, V+[x(t),t,dt] and V-\x{t),t, dt\. In order to go back to the 
“large-scale” domain and derive the “large-scale” velocities appearing in Eq. (2.16), we smooth out each frac¬ 
tal geodesic in the bundles selected by the zooming process with balls of radius larger than r. This amounts 
to carry out a transition from the non-differentiable to the differentiable domain and leads to define two “large- 
scale” velocity fields now resolution independent: V+[x{t),t,dt > r] = LS < V+[x{t),t,dt] >= v+[x{t),t] and 
V-[x{t),t,dt > r] = LS < V-[x{t),t,dt] >= v-[x{t),t]. The important new fact appearing here is that, after 
the transition, there is no more reason for these two velocities to be equal. While, in standard mechanics, the concept 
of velocity was one-valued, we must introduce, for the case of a non-differentiable space, two velocities instead of 
one, even when going back to the “large-scale” domain. This two-valuedness of the velocity vector finds its origin 
in a breaking of the discrete time reflection invariance symmetry (dt <-*■ —dt). Note that this symmetry is actually 
independent of the time reflection symmetry T [t ^ —t) of which the breaking will be considered, together with the 
breaking of the parity P (r ^ —r) and of the differential space symmetry (dr —dr), in Secs. tV and|^. 


Therefore, if one reverses the sign of the time differential element, becomes u_. But the choice of a reversed 
time element {—dt) must be as qualified as the initial choice (dt) for the description of the laws of nature, and we 
have no way, at this level of the description, to define what is the “ past” and what is the “future”. The only solution 
to this problem is to consider both the forward {dt > 0) and backward {dt < 0) processes on the same footing. Then 
the information needed to describe the system is doubled with respect to the standard differentiable description. 

A simple and natural way to account for this doubling consists in using complex numbers and the complex product. 
As we recall hereafter, this is the origin of the complex nature of the wave function of quantum mechanics, since 
this wave function can be identified with the exponential of the complex action that is naturally introduced in this 
framework. The choice of complex numbers to represent the two-valuedness of the velocity is not an arbitrary choice. 
The use of these numbers is the most natural way to generalize to two dimensions the set of real numbers. In the 
same way, as we shall see in Sec. IV, the use of bi-quaternionic numbers is the most natural way to generalize them 
to eight dimensions (see Appendix]^. 
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E. Covariant derivative operator 


Finally, we describe, in the scaling domain, the elementary displacements for both processes, dX±, as the sum of 
a large-scale part, dx± = z;± dt, and a fluctuation about this large-scale part, d^±, which is, by definition, of zero 
large-scale part, LS < d^± >= 0 

dX+{t) =v+dt + d^+{t), 

dX-(t) = V- dt + d^-{t). (2.25) 

Considering first the large-scale displacements, large-scale forward and backward derivatives, d/dt+ and d/dt-, are 
defined, using the large-scale part extraction procedure. Applied to the position vector, x, they yield the forward and 
backward large-scale velocities 


d , . d 

—x(i)=„+, —X(i)=»-. 


(2.26) 


As regards the fluctuations, the generalization to three dimensions of the fractal behavior of Eq. ( 2.20 ) writes (for 
D = 2) 


LS < d^±i d^±j >= ±2 T) Sij dt i, j = x, y, z, 


(2.27) 


as the d^(t)’s are of null large-scale part and mutually independent. The Kronecker symbol, Sij, in Eq. ( 2.27 ), implies 
indeed that the large-scale part of every crossed product d^±i d^±j, with i ^ j, is null. This is due to the fact 
that, even if each term in the product scales as dt^/"^, each behaves as an independent fractal fluctuation around its 
own large-scale part. Therefore, when we proceed to the smoothing out of the geodesic bundle during the transition 
from the small to the large-scale domain, we apply a process which is mathematically (not physically) equivalent to 
a stochastic “Wiener” process, and also more general, since we do not need any Gaussian distribution assumption. 
Thus, we can apply to the large-scale part of the d^{t) product the property of the product of two independent 
stochastic variables: i.e., the large-scale part of the product is the product of the large-scale parts, and therefore here 


zero. 


To recover local differential time reversibility in terms of a new complex process [[ll[ , we combine the forward and 
backward derivatives to obtain a complex derivative operator 


d 

dt 


d 

dt+ 


_d_ 

dt- 


d 

dt+ 


d 

dt- 


Applying this operator to the position vector yields a complex velocity 

d V++V- v+— V- 

V = —x(t) = V — lU = --- i --- 

dt ^ ' 2 2 


(2.28) 


(2.29) 


The minus sign in front of the imaginary term is chosen here in order to obtain the Schrodinger equation in terms 
of ijj- The reverse choice would give the Schrodinger equation for the complex conjugate of the wave function tjj'^, and 
would be therefore physically equivalent. 

The real part, V, of the complex velocity, V, represents the standard classical velocity, while its imaginary part, 
U, is a new quantity arising from non-differentiability. At the usual classical limit, Vj. = V- = v, so that V = v and 
U = 0. 


Contrary to what happens in the differentiable case, the total derivative with respect to time of a fractal function 
f{x{t),t) of integer fractal dimension contains finite terms up to higher order 


1 d^f dX,dXj 

dt dt dxi dt 2 dxidxj dt 


1 d^f dX.dXjdXk 
6 dxidxjdxk dt 


(2.30) 


Note that it has been shown by Kolwankar and Gangal [Q that, if the fractal dimension is not an integer, a 
fractional Taylor expansion can also be defined, using the local fractional derivative. 
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In our case, a finite contribution only proceeds from terms of H-order, while lesser-order terms yield an infinite 
contribution and higher-order ones are negligible. Therefore, in the special case of a fractal dimension D = 2, the 
total derivative writes 


^ 1 dXdX, 

dt dt dt 2 dxidxj dt 


Usually the t erm dXjdXj/dt is infinitesimal, but here its large- scale part reduces to LS < d^i d^j > /dt. Therefore, 
thanks to Eq. (p.27), the last term of the large-scale part of Eq. (p.31) amounts to a Laplacian, and we can write 


dt± 



- u±.V ± DA 


f ■ 


(2.32) 


Substituting Eqs. ( 2.32|) into Eq. ( ^.28 ), we finally obtain the expression for the complex time derivative operator 


d 

dt 


dt 


V.V - iVA . 


(2.33) 


The passage from standard classical (everywhere differentiable) mechanics to the new non-differentiable theory can 
now be implemented by a unique prescription: replace the standard time derivative d/dt by the new complex operator 
d/dt [^. In other words, this means that d/dt plays the role of a “covariant derivative operator” (in analogy with 
the covariant derivative DjA^ = djA^ + replacing djA^ in Einstein’s general relativity). 


F. Covariant mechanics induced by scale laws 


Let us now summarize the main steps by which one may generalize the standard classical mechanics using this 
covariance. We assume that the large-scale part of any mechanical system can be characterized by a Lagrange 
function C{x, V, t), from which an action S is defined 

pt2 

5 = / C{x,V,t)dt. (2.34) 

Jti 


Our Lagrange function and action are a priori complex and are obtained from the classical Lagrange function 
L{x, dx/dt, t) and classical action S precisely from applying the above prescription d/dt ^ d/dt. The action principle, 
which is a stationary action principle applied on this new complex action, leads to the generalized Euler-Lagrange 
equations Q 


d dC dC 
dt dVi dxi ’ 


(2.35) 


which are exactly the equations one would have obtained from applying the covariant derivative operator (d/dt 
d/dt) to the usual classical Euler-Lagrange equations themselves: this result demonstrates the self-consistency of 
the approach and vindicates the use of complex numbers (see Appendix ^). Other fundamental results of standard 
classical mechanics are also generalized in the same way. In particular, assuming homogeneity of space in the mean 
leads to define a generalized complex momentum and a complex energy given by 




(2.36) 


If we now consider the action as a functional of the upper limit of integration in Eq. (2.34), the variation of the 
action from a trajectory to another nearby trajectory, when combined with Eq. (2.35), yields a generalization of other 
well-known relations of standard mechanics 


7? = V5, E = -dS/dt. 


(2.37) 
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G. Generalized Newton-Schrodinger Equation 


Let us now specialize our study, and consider Newtonian mechanics, i.e., the general case when the structuring 
external field is a scalar potential, <&. The Lagrange function of a closed system, L = ^mv'^ — <&, is generalized, in the 
large-scale domain, as C{x, V, t) = The Euler-Lagrange equations keep the form of Newton’s fundamental 

equation of dynamics 

mi-V = -V$, (2.38) 

at 

which is now written in terms of complex variables and complex operators. 


In the case when there is no external field, the covariance is explicit, since Eq. ( 2.38 ) takes the form of the equation 
of inertial motion 


dV/dt^O. (2.39) 

This can be compared to what happens in general relativity, where the equivalence principle of gravitation and 
inertia leads to a strong covariance principle, expressed by the fact that one may always find a coordinate system 
in which the metric is locally Minkowskian. This means that, in this coordinate system, the covariant equation of 
motion of a free particle is that of inertial motion 


Du^ = 0, (2.40) 

where is the four-velocity of the particle and D, the covariant derivative operator of general relativity, defined by 
its application to any four-vector A'' as 


D^A'' = d^A'^ + T-^Ad 


Written in any coordinate system, Eq. (2.40|) becomes the local geodesic equation 


d^x>^ dx’^ dxP ^ 
ds^ ''P ds ds 


(2.41) 


(2.42) 


The covariance induced by scale effects leads to an analogous transformation of the equation of motions, which, as 
we show below, become the Schrodinger, Klein-Gordon and Dirac equations, and which we are therefore allowed to 
consider as local geodesic equations. 


In both cases, with or without external field, the complex momentum V reads 

V = mV, 


(2.43) 


so that, from Eq. (2.37), the complex velocity V appears as a gradient, namely the gradient of the complex action 

V = VS/m. (2.44) 

We now introduce a complex wave function '0 which is nothing but another expression for the complex action S 

= (2.45) 

The factor So has the dimension of an action (i.e., an angular momentum) and must be introduced at least for 
dimensional reasons. We show, in what follows, that, when this formalism is applied to microph ysics, So is nothing 
but the Ti fundamental constant. It is therefore related to the complex velocity appearing in Eq. (2.44) as follows 


V = -j —V(ln0). 


(2.46) 


We have now at our disposal all the mathematical tools needed to write the fundamental equation of dynamics of 
Eq. (2.38) in terms of the new quantity ip. It takes the form 


f5o^(Vln0) = V$. 


(2.47) 
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It is worth noting that the inertial mass m has disappeared from this equation. 


Now one should be aware that d and V do not commute. However, as we shall see in the followings, there is a 
particular choice of the arbitrary constant Sq for which d\ntjj)/dt is nevertheless a gradient. 


Replacing d/dt by its expression, given by Eq. ( 2.33|) , yields 

V$ = I5o ^^+V.V-zX>A^ (VlnV^), 


and replacing once again V by its expression in Eq. (p.46|), we obtain 


Vd) = iSvi 


d 


5, 


—V In V- - * <{ (V In V;.V)(V In V') + X>A(V In 


Consider now the remarkable identity 


(Vln/f + Aln/ = ^ 


which proceeds from the following tensorial derivation 
In / + 9^ In In / = 9^ ^ + 




f 

d^d^f 

f 


f f 


When we apply this identity to ip take its gradient, we obtain 


Alp 

iP 


= V[(Vln'(/i)^ + Ain'(/;]. 


(2.48) 


(2.49) 


(2.50) 


(2.51) 


(2.52) 


The second term in the right-hand side of this expression can be transformed, using the fact that V and A commute. 


i.e., 


VA = AV, 

and so can be the first term, thanks to another remarkable identity 

V(V/)2=2(V/.V)(V/), 


(2.53) 


(2.54) 


that we apply to / = liiip. We finally obtain 


V 



2(V In V)(V In ^/>) -k A(V In ^/>). 


(2.55) 


We recognize, in the right-hand side of this equation, the two terms of Eq. (^.49 ), 
of iSo/m and T>. Therefore, the particular choice 


which were respectively in factor 


5n = 2ml? 


(2.56) 


allows us to simplify the right-hand side of Eq. (2.49). The simplification is twofold: several complicated terms are 
compacted into a simple one and the final remaining term is a gradient, which means that the fundamental equation 
of dynamics can now be integrated in a universal way. The wave function in Eq. (2.45) is therefore defined as 


Ip = e 


iS j2m'D 


and is solution of the fundamental equation of dynamics, Eq. (2.38), which we write 

d . 


dt 


V = -222V 


+ 1 = - 


= —V^/m. 


(2.57) 


(2.58) 
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Integrating this equation finally yields 

V^/\'ib + iV^'ip-—il) = Q, (2.59) 

at 2m 

up to an arbitrary phase factor which may be set to zero by a suitable choice of the ip phase. 

We are now able to enlight the meaning of the choice Sq = 2m'D. At first sight, it is not a necessary condition 
from the viewpoint of physics, but merely a simplifying choice as regards the equation form. Indeed, it is only under 
this particular choice that the fundamental equation of dynamics can be integrated. However, if we do not make this 
choice, the tp function is a solution of a third order, non linear, complicated equation such that no precise physical 
meaning can be given to it. We therefore claim that our choice Sq = 2m'D has a profound physical significance, 
since, in our construction, the meaning of ip is directly related to the fact that it is a solution of the below obtained 
Schrodinger equation. 

In the case of standard quantum mechanics, as applied to microphysics, this result shows that there is a natural link 
between the Compton relation and the Schrodinger equation. In this case, indeed, Sq is nothing but the fundamental 
action constant h, while I? defines the fractal/non-fractal transition (i.e., the transition from explicit scale dependence 
to scale independence in the rest frame), A = 2I?/c. Therefore, the relation Sg = 2m'D becomes a relation between 
mass and the fractal to scale independence transition, which writes 

Ae = — . (2.60) 

me 

We recognize here the definition of the Compton lenght. Its profound meaning - i.e., up to the fundamental con¬ 
stants h and c, that of inertial mass itself - is thus given, in our framework by the transition scale from explicit 
scale dependence (at small scales) to scale independence (at large scales). We note that this lenght-scale should be 
understood as an “object” of scale space, not of standard space. 


We recover, in this case, the standard form of Schrodinger’s equation 


— Aip + ih—ip = <^ip. 


(2.61) 


The statistical meaning of the wave function (Born postulate) can now be deduced from the very construction of 
the theory. Even in the case of only one particle, the virtual geodesic family is infinite (this remains true even in the 
zero particle case, i.e., for the vacuum field). The particle properties are assimilated to those of a random subset of 
the geodesics in the family, and its probability to be found at a given position must be proportional to the d ensity 
of the geodesic fluid. This density can easily be calculated in our formalism, since the imaginary part of Eq. ( 2.59| ) 
writes in terms oi p = ipip^ 


dt 


-I- div(pl/) = 0, 


(2.62) 


where V is the real part of the complex velocity, which is identified, at the classical limit, with the classical velocity. 
This equation is recognized as the equation of continuity, implying that p = ipip^ represents the fluid density which is 
proportional to the probability density, thus ensuring the validity of Born’s postulate. The remarkable new feature 
that allows us to obtain such a result is that the continuity equation is not written as an additional a priori equation, 
but is now a part of our generalized equation of dynamics. 


III. COMPLEX KLEIN-GORDON EQUATION 


In Sec. H, the Schrodinger equation has been derived, from first principles, as a geodesic equation in a fractal three 
space, for non-relativistic motion in the framework of Galilean scale relativity. In the rest of this article we shall be 
concerned with relativistic motion in the same framework of Galilean scale relativity and shall be led to derive the 
corresponding free-particle quantum mechanical equations (Klein-Gordon’s and Dirac’s) as geodesic equations in a 
four-dimensional fractal space-time. 


The Klein-Gordon equation has already been established, in this framework, in a complex form ||l7| . We summarize 
and update, in the current section, the main steps of its derivation and shall be led to give a bi-quaternionic version 
of the same equation in Sec. IV. 
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A. Motion relativistic covariant derivative operator 


Most elements of the approach summarized in Sec. ^ remain correct in the motion-relativistic case, with the time, 
t, replaced by the proper time, s, as the curvilinear parameter along the geodesics. Now, not only space, but the full 
space-time continuum is considered to be non-differentiable, thus fractal. We consider a small increment dX^ of a 
non-differentiable four-coordinate along one of the geodesics of the fractal space-time. We can, as above, decompose 
dX^ in terms of a large-scale part LS < dX^ >= dx^ = Vfj,ds and a fluctuation respective to this large-scale part 
d^^, such that LS < d^^ >= 0, by definition. 

As in the non-relativistic motion case, the non-differentiable nature of space-time implies the breaking of the reflec¬ 
tion invariance at the infinitesimal level. If one reverses the sign of the proper time differential element, the large-scale 
part of the velocity v+ becomes V-. We therefore consider again both the forward {ds > 0) and backward {ds < 0) 
processes on the same footing. Then the information needed to describe the system is doubled with respect to the 
classical differentiable description. We have seen that this fundamental two-valuedness can be accounted for by the 
use of complex numbers. The new complex process, as a whole, recovers the fundamental property of microscopic 
reversibility. 


One is then led to write the elementary displacement along a geodesic of fractal dimension D = 2, respectively for 
the forward (-I-) and backward (—) processes, under the form 


dX^ = d±x^ + d^f^ = v^ds + u^V^ds^^'^, 


(3.1) 


with d±x^ = v'^ds and d^^ = . In these expressions, is a dimensionless fluctuation, and the length- 

scale 22? is introduced for dimensional purpose. We define the large-scale forw ard and backward derivatives, d/ds+ 
and d/ds-, using the large-scale part extraction procedure defined in Sec. |ll E|, as 


(3.2) 


*!/(“> =\-s;-/ ^ 


Once applied to x^, they yield the forward and backward large-scale four-velocities 


ds- 


-x'^{s) = 


ds- 


x^{s) = v^. 


(3.3) 


The forward and backward derivatives of Eq. (3.3) can be combined to construct a complex derivative operator 

(3.4) 


d 

ds 


d ^ d 
ds+ ds- 


d 

ds+ 


d 

ds- 


When applied to the position vector, this operator yields a complex four-velocity 

d vd + 

=—x>^ = V>^-iU>^ = ^ 
as 


2 2 

As regards the fluctuations, the generalization to four dimensions of the fractal behavior described in Eq. (| 

LS < d^^d^l >= T2Vr]i^'"ds. 


(3.5) 

3 ) gives 

(3.6) 


As noted in Sec. IIE , each term in the left-hand side product behaves as an independent fractal fluctuation around 
its own large-scale part, which justifies a mathematical treatement analogous to that of a stochastic Wiener process. 
We make in the present paper the choice of a (-I-, —, —, —) signature for the Minkoskian metric of the classical space- 
time 77 ^^. Therefore, the difficulty comes from the fact that a diffusion (Wiener) process makes sense only in 
where the “metric” 77 ^*^ should be positive definite, if one wants to interpret the continuity equation satisfied by the 
probability density as a Kolmogorov equation ||3^ . 

Several proposals have been made to solve this problem. Dohrn and Guerra introduce the above “Brownian 
metric” and a kinetic metric and obtain a compatibility condition between them which reads = 9^^- 

An equivalent method was developed by Zastawniak |^ , who introduces, in addition to the covariant forward and 
backward drifts and wO (be careful that our notations are different from his), new forward and backward drifts 
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and h^. In terms of these drifts the Fokker-Planck equations (that one can derive from the Klein-Gordon equation) 
become Kolmogorov equations for a standard Markov-Wiener diffusion process in R'^. Serva gives up Markov 
processes and considers a covariant process that belongs to a larger class, known as “Bernstein processes”. 

All these proposals are equivalent in the end, and amount to transform a Laplacian operator in into a Dalem- 
bertian. Namely, the two forward and backward differentials of a function f(x, s) write 

= (a 

As, in what follows, we only consider s-stationary functions, i.e., functions that do not explicitly depend on the 
proper time s, the complex covariant derivative operator reduces to 

^ + (3.8) 

The plus sign in front of the Dalembertian comes from the choice of the metric signature. 


B. Complex stationary action principle 


Let us now assume that the large-scale part of any mechanical system can be characterized by a complex action 
S. The same definition of the action as in standard relativistic motion mechanics allows us to write, according to the 
stationary action principle, 

SS = —me 5 j ds = 0, (3.9) 

J a 


between two fixed points a and 6, provided ds is defined as ds = LS < ^/dX^UiX^ >. Differentiating the integrand, 
we obtain 

5S = —me f Vy d{dx'^). (3.10) 

J a 


with 6x'' = LS < SX‘' >. 


Integrating by parts yields 


SS = -I 


Sx‘' 


ds 


(3.11) 


To get the equations of moti on, o ne has to determine <55 = 0 between the same two points, i.e., at the limits 
(Sx‘')a = (Sx‘')b = 0. From Eq. (3.11), we therefore obtain a differential geodesic equation 


dV, 

ds 


= 0 . 


(3.12) 


We can also write the elementary variation of the action as a functional of the coordinates. We have thus to consider 
the point a as fixed, so that {Sx'')a = 0. The second point b must be considered as variab le. T he only admissible 
solutions are those which satisfy the equations of motion. Therefore, the integral in Eq. (3.11) vanishes. Simply 
writing (Sx’')b as Sx'' gives 


SS = —me Vv Sx" 


(3.13) 


The minus sign in the right-hand side of Eq. ( 3.13 ) proceeds from our choice of the metric signature (see Sec. Ill A| ) . 
The complex four-momentum can thus be written as 


Vi, = meVv = —dvS. 


(3.14) 


Now, the complex action, 5, characterizes completely the dynamical state of the particle, and we can introduce a 
complex wave function 
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(3.15) 


_ JS/So 


tp = e 


It is linked to the complex four-velocity by Eq. (3.14), which gives 


Vjy = ^-^dylnip. 
me 


(3.16) 


C. Free-particle complex Klein-Gordon equation 


We now ap ply t he scale-relativistic prescription: replace the derivative in Eq. (3.12) by its covariant expression 
give n by Eq. (|3.8| ). Accounting for the expression of the complex four-velocity of Eq. ( 3.16 ) allows to transform 
Eq. (3.12) into 


C2 


■d^lnipdfidylnip — = 0 . 


(3.17) 


We see that the particular choice, Sq = Ti = 2mcV, anal ogous to the one made and discussed in Sec. IIG, allows 
us, once again, to simplify the left-hand side of Eq. (3.17), provided we make use of the following identity (which 
generalizes its three-dimensional counterpart given in Eq. (p.51)) 


2 V V' 


= ( d^lntp + -dfj, ) d^d’^lnip. 


(3.18) 


Dividing by the constant factor, we obtain the equation of motion of the free particle under the form 


d'' 




= 0 . 


(3.19) 


Therefore, the Klein-Gordon equation (without electromagnetic field) 

d^d^ip H-= 0, (3.20) 

h 

becomes an integral of motion of the free particle, provided the integration constant is chosen equal to a squared mass 
term, jf?. 

The quantum behavior described by this equation and the probabilistic interpretation given to tp is here reduced 
to the description of a free fall in a fractal space-time, in analogy with Einstein’s general relativity where a particle 
subjected to the effect of gravitation is described as being in free fall in a curved space-time. 


Moreover, these equations can indeed be considered as “covariant”, since the relativistic quantum equation written 
in terms of the complex derivative d/ds has the same form as the equation of a relativistic macroscopic and free 
particle written in terms of the usual derivative d/ds. 


It is worth noting here that, contrary to the hope expressed in Ref. |28|, the metric form of special and general 
relativity, = I, is not conserved in quantum mechanics. Indeed, it has been shown by Pissondes that the 

free particle Klein-Gordon equation, when it is expressed in terms of the complex four-velocity V^, leads to the new 
following equality 


-f = 1. 


(3.21) 


In the scale-relativistic framework, this expression defines the metric that is induced by the internal scale structures 
of the fractal space-time. 

In the absence of an electromagnetic field, and S are related by Eq. (|3.14), which we can write 


V. =- d.S, 


me 


(3.22) 


17 






















so that Eq. (3.21) becomes 


d'^Sdf_iS — 2imc'Dd^d^S = rn^c^, 


(3.23) 


which is the new form taken by the Hamilton-Jacobi equation. 


Another important property, which comes to light in the above derivation but was hidden in Ref. |17|, is the way 
the characteristic length of the fractal (quantum) to differentiable (classical) transition (the Compton length of the 
particle in its rest frame) and the mass term in the Klein-Gordon equation appear. Here, as in the non-relativistic 
motion case (Schrodinger), we identify the transition scale 27) to the Compton length 7i/mc, provided Sq is the 
fundamen tal ac tion constant Ti, in order to obtain the motion equation under the form of a vanishing four-gradient, 
as in Eq. (|3.19|) . Should we make another choice, we would not be able to integrate this equation and would obtain a 
third order, non linear, complicated one to which no precise physical meaning could be given. As for the mass term 
in the final Klein-Gordon equation, it does not appear to proceed from a choice of the transition length, as it was 
claimed in Ref. [p^, but to correspond to a mere integration constant. 


IV. BI-QUATERNIONIC KLEIN-GORDON EQUATION 

It has long been known that the Dirac equation naturally proceeds from the Klein-Gordon equation when written 
in a quaternionic form [}l2|j43[ . We propose in the current section to introduce naturally a bi-quaternionic covariant 
derivative operator, leading to the definition of a bi-quaternionic velocity and wave-function, which we use to derive 
the Klein-Gordon equation in a bi-quaternionic form. This allows us to obtain, in Sec.^, the Dirac equation as a mere 
consequence, since Dirac spinors and bi-quaternions are actually two equivalent representations of an electron. 

The quaternionic formalism, as introduced by Hamilton |4^, and further developed by Gonway |43y45| ] is recalled 
in Appendix p. 


A. Symmetry breaking and bi-quaternionic covariant derivative operator 


Most of the approach described in Sec. HI remains applicable. However, the main new features obtained in the 
now studied case proceed from a deeper description of the scale formalism, considering the more general case when 
the peculiar choice of an identification of the differentials and the resolution variables (remember the choice explicited 
in Sec. H C, where we have set dt = 6t, which corresponds, in the motion relativistic case, to ds = Ss) is given up, 
implying the subsequent breaking of the symmetries: 
ds —ds 
dx>^ ^ —dx^ 


We have already stressed, in Sec. H.C, that differentiability can be recovered (at large scales), even when dealing 
with non-differentiability. However, this implies that, for any set of differential equations describing a given process, 
the physical and mathematical descriptions are only coincident on a limited scale range. We can say that any con¬ 
sistent mathematical tool lives in a description space which is tangent to the physical space and that the validity 
of this tool is therefore limited to a finite scale region around the contact point. For the Schrodinger and complex 
Klein-Gordon equations the tangent mathematical space is such that dt = 6t (Schrodinger) or ds = Ss (Klein-Gordon). 
When jumping to smaller (higher energy) scales, we are led to give up this peculiar choice, and retain, at the Dirac 
scale, a new mathematical description (detailed below) where the differentials and the resolution variables do no more 
coincide. However, we are aware that this description is only valid at the scales where the Dirac equation applies 
without corrections and that we shall be led to improve it when going to yet smaller scales. We know actually that, 
at higher energy, the probability amplitudes must become more complicated to take into account, e.g., the isotopic 
spin, the hypercharge, etc.; moreover, one must take into account the radiative corrections which, in scale relativity, 
are identified to new scale dependent internal structures, through their description in terms of the renormalisation 
group. This scale depending nature of the coincidences between the mathematical and physical worlds can be viewed 
as a truism since it is implicit in the ordinary run of physics, but it becomes a key character of the scale-relativistic 
formalism where the choice of the differential variables must be explicited at each stage of the theory. 
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In the scaling domain, the four space-time coordinates X^(s, e^, e^) are four fractal functions of the proper time s 
and of the resolutions for the coordinate and for the proper time. Q We consider the case when, for an elemen¬ 
tary displacement corresponding to a shift ds in the curvilinear parameter, the resolutions verify < dX^ and 
Eg < ds, which implies that, at a given a forward shift ds of s yields a displacement dX^ of X^ and a backward 
shift —ds produces a displacement —dX^, of which the amplitudes are not necessarily equal (see Fig. ^). 



We can therefore apply, to these two different elementary displacements, the canonical decomposition 

dXf^ = dx>^+d^>^, (4.1) 

with 

TS < dX>^ >= dx^ = v^+ ds, (4.2) 


and 


with 


dC = alV^ids^)^, LS <al>=0, LS < {a'l f >= 1, 


- dX>^ = Sx^ + 6^^, 


LS < -dX>^ >= 5x>^ = vt ds. 


(4.3) 

(4.4) 

(4.5) 


= atV^{ds^)^, LS<at>=Q, LS<{a^_f>=l. (4.6) 

In the differentiable case, dX^^ = —{—dX^), and therefore = —v^_. This is no longer the case in the non- 
differentiable case, where the local symmetry dx^ —dx^ is broken. 

Furthermore, we must also consider the breaking of the symmetry ds —ds proceeding from the twofold definition 
of the derivative with respect to the curvilinear parameter s. Applied to A^, considering an elementary displacement 


^Please note that the e^’s are not actnally vectors but the diagonalized version of a tensorial quantity e^i,, formally analogous 
to an error ellipsoid. This point, irrelevant to the present study, will be developed in forthcoming works. 
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dX^, the large-scale part extraction process gives two large-scale forward and backward derivatives d/ds+ and d/ds-, 
which yield in turn two large-scale velocities, which we denote +. Considering now the same extraction process, 

s 

applied to an elementary displacement —dX^^ the large-scale forward and backward derivatives d/ds+ and d/ds- 
allow again to define two large-scale velocities, denoted v± _ . We summarize this result as 

s /i, 



dx^^ 
ds± ’ 


V± - 


s fj. 


Sx^ 

ds± 


(4.7) 


In Galilean scale relativity, with a constant fractal dimension D = 2, we can, at this stage, define several total 
derivatives with respect to s of a fractal function /. We write them - using a compact strai ghtfo rward notation with 
summation over repeated indices - after substituting, in the four-dimensional analog of Eq. (2.31), the expressions for 
the derivatives of the X'^ obtained when using Eqs. (4.1) to ( ^(^ ) with the different expressions of Eq. (4.7) for the 
large-scale velocities 


df « u N 5/ u dFf 

ds ‘ ^ V ^ t as a ij. a ^ oX^^ dXt^dX'^ 


with 


w 


= a'"y^(ds^) 2D 2. 


(4.8) 


(4.9) 




Now, when we apply the large-scale part operator to Eq. ([4.8]), using Eq. ( pd)| ) and the properties of the a5_ 's as 
stated in Eqs. (O) and (4^), the ’s disappear at the first order, but, at the second order, for the fractal dimension 
U = 2, we obtain the analogue of Eq. (^) 


LS < ±w'± ± >= ^2'Dt]^'^ ds, 

S s fj, 

the =F sign in the right-hand side being the inverse of the s-sign in the left-hand side. 


(4.10) 


We can therefore write the large-scale part of the total derivatives of Eq. (4.8) as 


ds ^ V ^ t \as a IJ. J 

where the =F sign in the right-hand side is still the inverse of the s-sign. 

When we apply these derivatives to the position vector X^, we obtain, as expected, 

dX^^ _ ^ 

^ ^ ± • 

as s M s u 


(4.11) 


(4.12) 


We consider now the four fractal functions —X^(s, e^, Cg). At this description level, there is no reason for 
(—W^)(s, e^, Cg) to be everywhere equal to —(X^)(s, e^, Cg), owing to a local breaking of the P (for fi = x,y,z) 
and T (for y = t) symmetries. Furthermore, as we have stressed above for X^, at a given —X^, a forward shift ds of 
s yields a displacement d{—X^) of —X^ and a backward shift —ds produces a displacement —d{—X^), with no nec¬ 
essarily equal amplitudes. We can therefore apply, to these two different elementary displacements, a decomposition 
similar to the one described in Eqs. ([4.lD to (4.6), i.e., 


with 


d{-Xf^) =dx>^ +de, 

TS < d{-Xf^) >= ds, 


(4.13) 

(4.14) 


d^>^ = , LS <d>l>=0, 


LS < (d^y >= 1, 


(4.15) 


and 
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d{-Xf^) = 6x^ + 


(4.16) 


with 


LS < -d(-X^) >= = vt ds 


(4.17) 


= d^V^idt^ 


LS < at >= 0, LS < {dty >= 1, 


(4.18) 


Then, we jump to the large-scale parts and are once more confronted to the breaking of the ds —ds symmetry. 
Applied to —X^, considering an elementary displacement d{—X^), the large-scale forward and backward derivatives 
d/ds+ and d/ds- give again two large-scale velocities, which we denote vt + ■ After another large-scale part extraction, 

considering now an elementary displacement —d(—X^), the large-scale forward and backward derivatives d/ds+ and 
d/ds- yield two other large-scale velocities, denoted v± - ■ We write in short 


+ 


dx>^ 
ds± ' 




Sx^ 

ds± 


(4.19) 


We therefore obtain new different total derivatives with respect to s of a fractal function /, which we write, as 
above, 


df Of 

— ±±±±± = -—-|-(u 

ds ^ ^ y ^ i OS s tj. 


M d^f 

±± + w± ± )-H-^ + at-a^V - 


^X^^ 


dX^^dX'^ 


with 


w‘ 


= d>^Vw{dsy^-X 


(4.20) 


(4.21) 


The large-scale part operator applied to Eq. ( 4.20| ), as done in Eq. (4.11) for Eq. (4.8), yields the large-scale total 
derivatives 


=(^+it±d^TVd>^d^]f. 

ds o X B ^ t \(js s IJ. 


(4.22) 


When we apply these derivatives to the position vector A^, we obtain, as expected again, 


ds 


■± ± = ± . 


(4.23) 


If we consider the simplest peculiar case when the breaking of the symmetry dx^ —dx^ is isotropic as regards the 
four space-time coordinates (i.e., the signs corresponding to the four /i indices are chosen equal), we are left with eight 
non-degenerate components - four vt ± and four ± - which can be used to define a bi-quaternionic four-velocity. 

s fj, s 

We write 


V'^=7r( 




0_)--«+-D(l_) + 

\ivt_ + v>l_) - - vt+) 




ei 


62 








63- 


(4.24) 


The freedom in the choice of the actual expression for will be discussed l ater. It is constrained fy the following 
requirements: at the limit when —> dX^ and Cg —*■ ds, every ei-term in Eq. (4.24) goes to zero, and, as = vt^ 
in this limit, one recovers the complex velocity of Eq. (2.29), -I- vt+ — *(i'++ — vt+)]/2-, at the classical 

"++■ 


limit, every term in this equation vanishes, save the real term, and the velocity becomes classical, i.e., real: = vt 


The bi-quaternionic velocity thus defined corresponds to a bi-quaternionic derivative operator (f/ds, similarly de¬ 
fined, and yielding, when applied to the posit ion vector A^, the corresponding velocity. Eor instance, the derivative 
operator attached to the velocity in Eq. (4.24) writes 
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++ 


ds 



+- 


-+ 


ds 


-+ 


o -r+- - -r++ ei 


ds 


++ 



ds 


-+ 


ds 


+- 


63- 


(4.25) 


Substituting Eqs. (4.11) and (4.22) into Eq. (4.25), we obtain the expression for the bi-quaternionic proper-time 
derivative operator 


d 


d 


— = [1 -I- ei -I- 62 -I- (1 - ^ 63 ]— -I- + iVd^di 


ds - ' - ■ — - —(4-26) 

the -I- sign in front of the Dalambertian proceeding from the choice of the metric signature (-1-,—,—,—). We keep 
here, for generalty, the d/ds term, stressing that it is actually of no use, since the various physical functions are not 
explicitly depending on s. It is ea sy to check that this operator, applied to the position vector X^, gives back the 
bi-quaternionic velocity of Eq. ([1.24 ). 


It is worth noting that the expression we have written for in Eq. (4.24) is one among the various choices we 
could have retained to define the bi-quaternionic velocity. The main constraint limiting this choice is the recovery 
of the complex and real velocities at the non-relativistic motion and classical limits. We also choose such as 


to obtain the third term in the right-hand side of Eq. (4.26) under the form of a purely imaginary Dalembertian, 


which allows to recover an integrable equation of motion. To any bi-quaternionic velocity satisfying both prescriptions 
corresponds a bi-quaternionic derivative operator (f/ds, similarly defined, and yielding this velocity when applied to 
the position vector X^. But, whatever the definitio n reta ined, the derivative operator keeps the same form in terms 
of the bi-quaternionic velocity V^, as given by Eq. ([4.32 ). Therefore, the different choices allowed for its definition 
merely correspond to different mathematical representations leading to the same physical result. 


B. Bi-quaternionic stationary action principle 


We now apply the stationary action principle, as developed in Sec. IIIB, to an elementary var iatio n of a bi- 

(4.27) 


quaternionic action S. The free motion equation issued from this principle has the same form as Eq. (3.12), i.e., 

dV. 


ds 


= 0 , 


where is the bi-quaternionic four-velocity, e.g., as defined in Eq. (4.24). 

The elementary variation of the action, considered as a functional of the coordinates, gives once more 

SS = —me Vfj, Sx^. 

We thus obtain the bi-quaternionic four-momentum, as 

= meV^ = -d^S. 

We can now introduce a bi-quaternionic wave function, which is again a re-expression of the action and which we 
write 


(4.28) 


(4.29) 


'-dull) = 


daS, 


cSo ^ 

using, in the left-hand side, the quaternionic product defined in Eq. 
velocity, as derived from Eq. (4.29), 

m 


(4.30) 

|). This gives for the bi-quaternionic four- 

(4.31) 


It is worth stressing here that we could choose, for the definition of the wave function in Eq. ( 4.3C ), a commutated 
expression in the left-hand side, i.e., {dfj,'ip)'4’~^ instead of ip~'^dfj,ip. But with this reversed choice, owing to the non- 
commutati vity of the quaternionic product, we could not obtain the motion equation as a vanishing four-gradient, 
as in Eq. (4.38). Therefore, we retain the above choice as the simplest one, i.e., yielding an equation which can be 
integrated. 
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C. Free-particle Klein-Gordon equation 


Our next step will be to apply the scale-relativistic presc ripti on: replace the standard proper time derivative d/ds 
by the new bi-quaternionic operator (f/ds, defined in Eq. (4.26). As, in what follows, we only consider s-stationary 
functions, i.e., functions which do not explicitly depend on the proper time s, the derivative operator reduces to 


(f 

ds 


= V''d,+iVd'^d,. 


(4.32) 


Now this expression is independent of the peculiar choice retained for t he bi-quaternionic form of the four-velocity, 
as the representation dependent term in the right-hand side of Eq. (4.26|) has vanished. 


The equation of motion, Eq. (4.27), thus writes 

(v"a, + ivd^'d,) = 0. 


We replace respectively V*", by their expressions given in Eq. ([4.3l|) and obtain 

i— (d''ijjdv + i'Dd''di\ = 0. 

m \ m J 


As in Sec. 11G, the choice 5o = 2mV allows us to simplify this equation and get 




The definition of the inverse of a quaternion 

= 1 , 


(4.33) 


(4.34) 


(4.35) 


(4.36) 


implies that tjj and commute. But this is not necessarily the case for ip and nor for iIj~^ and d^ip and their 

contravariant counterparts (see Appendix ^). However, when we derive Eq. (4.36) with respect to the coordinates, 
we obtain 


Ip df^ip ^ = -{d^ip)ip ^ 

ip~'^d^ip = -{d^ip~^)ip, (4.37) 

and identical formulae for the contravariant analogues. 


Developing Eq. (4.35), using Eqs. (4.37) and the property d'^d^d^ = d^d^d^, we obtain, after some calculations, 


df,[{d’'d,iP)ip-^] = 0. (4.38) 

We integrate this four-gradient and write 

{d''d,iP)i>-^ +C = 0, (4.39) 

of which we take the right product by ip to obtain 

d'^d^ip + Cip = 0. (4.40) 

We therefore recognize the Klein-Gordon equation for a free particle with a mass m so that rr?(? = C. 


V. DIRAC EQUATION 

We now use a long-known property of the quaternionic formalism, which allows to obtain the Dirac equation for a 
free particle as a mere square root of the Klein-Gordon operator (see, e.g.. Refs. |^,Q). We give below a pedestrian, 
but accurate, derivation of this property. 

We first develop the Klein-Gordon equation, as 
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(5.1) 




c^ dt^ dx^ dy"^ 


dz"^ 


Thanks to the property of the quaternionic and complex imaginary units = e| = e| = = — 1, we can write 

Eq. (5.1) under the form 


1 d'^ip 2^^ 2 , • 219^ . , 2^^V' 2 , -2 

+ i 


^ 2/2 
^2 


(5.2) 


We now take advantage of t he anticommutative property of the quaternionic units {eiSj = —ejSi for i ^ j) to add 
to the right-hand side of Eq. (pj) six vanishing couples of terms which we rearrange to obtain 


Id fl dll: 


cdt \c dt 


9 / 

dip 

dip . 


e3T;-e2 + ei 
^ dx 


9 / 

dip 

90. 

'dy ( 

63 - 5-62 + 61 
^ OX 

9j/* 

9 / 

dip 

90, 

^^Yz[ 

63T;-62 + 61 

ox 

dy^ 

me 

( dih 

dip 

Ti 

ea — 62 -hei — 
\ ox oy 


dip me \ 
ea-^ei - i — e^ipe^ 1 62 

dip me \ 
63-^61 - i — esipes j i 


dip me 

iTr^i - *-^63 
dz n 

dip me 

S 3 ^ei - i—t 

dz n 


(5.3) 


We see that Eq. (5.3) is obtained by applying twice to the bi-quaternionic wavefunction ip the operator d/edt 
written as 


19 d d d me , , 

“"dl = ®3^e2-b ei—i-b 63 —ei - i—63) )e3. 

c dt dx dy dz n 


(5.4) 


The thre e first Conway matrices e 3 ( )e 2 , ei( )i and e 3 ( )ei (see Appendix ^), figuring in the right-hand side 

of Eq. (5^), can be written in the compact form — with 


= 


0 (Tfe 
(Tfe 0 


the (Tfe’s being the three Pauli matrices, while the fourth Conway matrix 

eai )e3 = 



can be recognized as the Dirac j3 matrix. We can therefore write Eq. (^) as the non-covariant Dirac equation for a 
free fermion 


1 9-0 k dip . me ^ , 

The covariant form, in the Dirac representation, can be recovered by applying 163 ( )e 3 to Eq. (^.5|). 


(5.5) 


The isomorphism which can be established between the quaternionic and spinorial algebrae through the multipli¬ 
cation rules applying to the Pauli spin matrices allows us to identify the wave function 0 to a Dirac spinor. Spinors 
and quaternions are both a representation of the SL(2,C) group. See Ref. Q for a detailed discussion of the spinorial 
properties of bi-quaternions. 


VI. CONCLUSION 

Three fundamental motion equations which are merely postulated in standard microphysics have been recovered, 
in the framework of Galilean scale relativity, as geodesic equations in a fractal space (Schrodinger), then in a frac¬ 
tal space-time (free Klein-Gordon and Dirac). It is interesting to note how the change from classical to quantum 
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non-relativistic motion, then from quantum non-relativistic to quantum relativistic motion arises from successive 
symmetry breakings in the fractal geodesic picture. 


First, the complex nature of the wave function is the result of the differential (proper) time symmetry breaking, 
which is the simplest effect arising from the fractal structure of space (space-time). It allows a statistical interpre¬ 
tation of this wave function in the form of a Born postulate. This interpretation is here no more a postulate, since 
the continuity equation we write in above Sec. 11G is now a part of our derived equation of dynamics. At this stage, 
Galilean scale relativity with a complex wave function permits the recovery of both the Schrodinger (Sec. ||) and 
Klein-Gordon (Sec. Ill) equations. 


To go on with the description of the elementary properties encountered in the microphysical world, we have to 
consider further breakings of fundamental symmetries, namely the differential coordinate symmetry {dx^ <-> —dx^) 
breaking and the parity and time reversal symmetry breakings. These new breakings provide a four-complex compo¬ 
nent wave function (i.e., a eight component wave function), of which the most natural mathematical representation is 
in term of bi-quaternionic numbers (see Appendix |a|). We therefore obtain the spinorial and the particle anti-particle 
nature of elementary objects which we can describe as Dirac spinors. Here, spin arises from the isomorphism between 
the quaternionic and spinorial representations, both of which are different representations of the SL(2,C) group. At 
this stage, the Klein-Gordon equation can be recovered from a mere stationary action principle in a bi-quaternionic 
form which naturally yields the free Dirac equation. 


It is worth stressing that these important results only proceed from a very restricted use of the scale-relativistic 
potentialities. We have, in the present work, limited our investigations to the induced effects of scale laws on the 
equations of motion, in the framework of dilation laws exhibiting a Galilean group structure, i.e., a fractal space-time 
with a constant fractal dimension D = 2. As it has been explained in the introduction, this is only one of the simplest 
levels at which the scale-relativistic program can be realised. We have also made a series of simplifying choices, which 
we have explained and justihed all along the derivation procedure. Some of them have been dictated by physical or 
experimental considerations, but others only correspond to special cases, provisionally retained, and other possibilities 
will have to be explored in the future. We keep therefore confidence in the fact that the further developments of the 
theory should be able to bring new insights into the understanding of the laws of nature. 


APPENDIX A: COMPLEX NUMBERS AND BI-QUATERNIONS AS A COVARIANT 

REPRESENTATION 


1. Introduction 


If one attempts to trace back the origin of most of the quantum behavior, one finds it in the complex nature of the 
probability amplitude. Feynman particularly stressed this point |^ , noticing that, once the first axiom of quantum 
mechanics is admitted (i.e., the probability amplitude ip is complex, it is calculated as tp = tpi + tjj 2 for two alternative 
channels, and the probability density is given by ipip^), quantum laws as well as classical laws can be naturally re¬ 
covered. Is it possible to “understand” this axiom? In other words, is it possible to derive it from more fundamental 
principles, for which an intuitive comprehension would be possible? 


It is clear that such an understanding is impossible in the framework of quantum mechanics itself (since these 
statements are its basic axioms), but must be looked for in an enlarged paradigm. The theory of scale relativity 
provides us with an extension, both of the foundation of physical theories and of the principle of relativity. So it could 
be worth asking such questions in its frame of thought. In scale relativity, we extend the founding stones of physics 
by giving up the hypothesis of space-time differentiability; then we extend the principle of relativity by applying it, 
not only to motion laws, but also to scale laws. 


We have recalled in Sec. IID our demonstration 0 that one of the main and simplest consequences of the non¬ 
differentiability of space-time is that the velocity vector becomes two-valued. This two-valuedness of the velocity 
implies a two-valuedness of the Lagrange function, and therefore a two-valuedness of the action S. Finally, the wave 
function is defined as a re-expression for the action, so that it will also be two-valued in the simplest case. 
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But we have, up to now, admitted without justification that this two-valuedness is to be described in terms of 
complex numbers. Why? This is an important point for the understanding of quantum mechanics, but also from the 
viewpoint of the meaning of the scale-relativistic approach. Indeed, our equations are not simply a “pasting” of real 
and imaginary equations, but involve the complex product from the very beginning of the calculations. In particular, 
the geodesic equation for a free particle takes the form of the equation of inertial motion 



= 0 , 


(Al) 


in terms of the “covariant” complex time derivative operator 


d 

dt 


^ + V.V -iVA. 
at 


(A2) 


Equation ( |A1| ) amounts to Schrodinger’s equation (see main text). Since it corresponds to a second derivative, the 
complex product has mixed the real and imaginary quantities in a very specific way, so that it is not at all a trivial 
statement that Schrodinger’s equation is recovered in the end. 


The aim of this appendix is precisely to address this specific question: why complex numbers and why bi¬ 
quaternions? As we shall see, the answer is that complex numbers achieve a particular representation of quantum 
mechanics in terms of which the fundamental equations take their simplest form. Other choices for the represen¬ 
tation of the two-valuedness and for the new product are possible, but these choices would give to the Schrodinger 
equation a more complicated form, involving additional terms (although its physical meaning would be unchanged). 
Bi-quaternions follow as a further splitting, at another level, of the complex numbers. 


2. Origin of complex numbers and quaternions in quantum mechanics 

Let us return to the step of our demonstration where complex numbers are introduced. All we know is that each 
component of the velocity now takes two values instead of one. This means that each component of the velocity 
becomes a vector in a two-dimensional space, or, in other words, that the velocity becomes a two-index tensor. So let 
us introduce generalized velocities 


V^ = {V\U^) ■ fc= 1,2,3 ; a = -,+. (A3) 

This can be generalized to all other physical quantities affected by the two-valuedness: namely, scalars A of the 
position space become vectors of the new 2D-space, etc... . While the generalization of the sum of these quantities 
is straighforward, = A^ -|- the generalisation of the product is an open question, which should be derived from 
first principles. 


The problem can be put in a general way: it amounts to find a generalization of the standard product that keeps 
its fundamental physical properties, or at least that keeps the most important ones (e.g., internal composition law) 
when some of them cannot escape to be lost (e.g., commutativity when jumping to quaternions). 


From the mathematical point of view, we are here exactly confronted to the well-known problem of the doubling of 
algebra (see, e.g., Ref. Q). Indeed, the effect of the symmetry breaking dt —dt (or ds <-*■ —ds) is to replace the 
algebra A in which the classical physical quantities are defined, by a direct sum of two exemplaries of A, i.e., the space 
of the pairs (a, b) where a and b belong to A. The new vectorial space A^ must be supplied with a product in order to 
become itself an algebra (of doubled dimension). The same problem is asked again when one takes also into account 
the symmetry breakings dx —dx and x^ —x^, which is the aim of the present paper: this leads to new algebra 
doublings. The mathematical solution to this problem is well-known: the standard algebra doubling amounts to 
supply A^ with the complex product. Then the doubling IR^ of IR is the algebra® of complex numbers, the doubling 
of® is the algebra IH of quaternions, the doubling of quaternions is the algebra of Graves-Cayley octonions. 
The problem with algebra doubling is that the iterative doubling leads to a progressive deterioration of the algebraic 
properties. Namely, the quaternion algebra is non-commutative, while the octonion algebra is also non-associative. 
But an important positive result for physical applications is that the doubling of a metric algebra is a metric algebra 
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These mathematical theorems fully justify the use of complex numbers, then of quaternions, in order to describe 
the successive doublings due to discrete symmetry breakings at the infinitesimal level, which are themselves more and 
more profound consequences of space-time non-differentiability. 

However, we give in what follows complementary arguments of a physical nature, which show that the use of the 
complex product in the first algebra doubling have a simplifying and covariant effect (we use here the word “covariant” 
in the original meaning given to it by Einstein namely, the requirement of the form invariance of fundamental 
equations). 

In order to simplify the argument, let us consider the generalization of scalar quantities, for which the product law 
is the standard product in IR. 

The first constraint is that the new product must remain an internal composition law. We also make the simplifying 
assumption that it remains linear in terms of each of the components of the two quantities to be multiplied. Therefore 
we consider a general form for a bilinear internal product 

(A4) 

where the matrix is a tensor (similar to the structure constants of a Lie group) that dehnes completely the new 
product. 


The second physical constraint is that we recover the classical variables and the classical product at the classical 
limit. The mathematical equivalent of this constraint is the requirement that A still be a sub-algebra of A?. Therefore 
we identify ao € A with (ao, 0) and we set (0,1) = a. This allows us to write the new two-dimensional vectors in the 
simplihed form a = qq + aia, so that the product now writes 

c = (ao + aia) (bo + bia) = aobo -I- aibia^ + (aobi + aibo)a. (AS) 

The problem is now reduced to find a^, i.e., only two fl coefficients instead of eight 

= ujo + ujia. (A6) 


Let us now come back to the beginning of our construction. We have introduced two elementary displacements, 
forward and a backward one, each of them made of two terms, a large-scale part and a fluctuation (see Eq. (2.25|)) 


dX+(t) = v+ dt + d^+(t), 

dX-(t) = V- dt + d^-(t). (A7) 


Therefore, one can define velocity fluctuations = d^^/dt and = d^-/dt, then a complete velocity in the 
doubled algebra 1^ 


V + W = 


f v+ + V- 

V 2 



f w+ -I- W- 

V 2 



(A8) 


Note that in terms of standard methods, this writing would be forbidden since the velocity W is infinite from 
the viewpoint of usual differential calculus (it is oc dt~^^'^). But we recall that we give meaning to this concept by 
considering it as an explicit function of the differential element dt, which becomes itself a variable. 


Now, from the covariance principle, the Lagrange function in the Newtonian case should strictly be written: 

C = ^mLS < (V + Wf >= {TS <V^ > +LS < >) (A9) 

We have LS < W >= 0, by definition, and LS < VW >= 0, because they are mutually independent. But what 
about LS < yV^ > 1 The presence of this term would greatly complicate all the subsequent developments toward 
the Schrodinger equation, since it would imply a fundamental divergence of non-relativistic quantum mechanics (note 
however that such a divergence finally happens in relativistic quantum held theories and leads to the renormalisation 
and renormalisation group approachs). Let us develop it 

ALS < yV^ > = LS < [(w+ + W-) — a (w+ — rc-)]^ > 

= LS < (w^ -I- ra?.)(l -I- a^) — 2a(w^ — vA_) + 2w+w-(l — a^) > . (AIO) 
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Since LS < vo]^ >= LS < vo^ > and LS < w+w- >= 0 (they are mutually independent), we finally find that 
LS < > can only vanish provided 


a^=-l, (All) 

namely, a = i, the imaginary. Therefore we have shown that the choice of the complex product in the algebra 
doubling plays an essential physical role, since it allows to suppress what would be additional infinite terms in the 
final equations of motion. 


3. Origin of bi-quaternions 


A last point that must be justified is the use of complex quaternions (bi-quaternions) for describing the new two- 
valuedness that leads to bi-spinors and the Dirac equation. One could think that the argument given in Sec. A 2 
(algebra doubling) leads to use Graves-Cayley octonions (and therefore to give up associativity) in the case of three 
successive doublings as considered in this paper. However, these three doublings are not on the same footing from a 
physical point of view: 

(i) The first two-valuedness comes from a discrete symmetry breaking at the level of the differential invariant, namely, 
dt in the case of a fractal space (yielding the Schrodinger equation) and ds in the case of a fractal space-time (yielding 
the Klein-Gordon equation). This means that it has an effect on the total derivatives d/dt and d/ds. This two- 
valuedness is achieved by the introduction of complex variables. 

(ii) The second two-valuedness (differential parity and time reversal violation, “dX”), which is specifically introduced 
and studied in the present paper, comes from a new discrete symmetry breaking (expected from the giving up of 
the differentiability hypothesis) on the space-time differential element dx^ <-^ —dx^. It is subsequent to the first 
two-valuedness, since it has an effect on the partial derivative d/j, = d/dx^ that intervenes in the complex covariant 
derivative operator, namely, 


^ + iVd^^)d^. (A12) 

ds 

(iii) The third two-valuedness is a standard effect of parity (“P”) and time reversal (“T”) in the motion-relativistic 
situation, which is not specific of our approach and is already used in the standard construction of Dirac spinors. 
It does not lead to a real information doubling, since Dirac spinors still have only four degrees of freedom as Pauli 
spinors do. 


Therefore, from the second and third doublings, complex numbers, then quaternions, can be introduced, which 
will affect variables which are already complex due to the first, more fundamental, doubling. This leads to the bi- 
quaternionic tool we use in the present work. 


Let us conclude this section by noting that these symmetry breakings are effective only at the level of the underlying 
description of elementary displacements (namely, in the non-differentiable fractal space-time). The effect of introducing 
a two-valuedness of variables in terms of double symmetrical processes precisely amounts to recover symmetry in terms 
of the bi-process, and therefore in terms of the quantum tools which are built from it. But reversely, this remark opens 
a possible way of investigation for future research about the origin of other features characteristic of the microphysical 
world. 


4. Complex representation and covariant Euler-Lagrange equations. 

Let us confirm by another argument that our choice of a representation in terms of complex numbers for the two- 
valuedness of the velocity (namely, in the minimal situation that leads to the Schrodinger equation) is a simplifying 
and covariant choice. We indeed demonstrate hereafter that the standard form of the Euler-Lagrange equations is 
conserved when we combine the forward and backward velocities in terms of a unique complex velocity. 

In a general way, the Lagrange function is expected to be a function of the variables x and their time derivatives 
X. We have found that the number of velocity components x is doubled, so that we are led to write 

L = L{x,x+,x-,t). (A13) 
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Instead, we have made the choice to write the Lagrange function as 


L = L{x,V,t). 


(A14) 


We now justify this choice by the covariance principle. Re-expressed in terms of x+ and the Lagrange function 
writes 


Therefore we obtain 


rl 1+*. , 

L = L[x, x+ + x-,t] . 


dL 1 — i dL dL 1 + i dL 


dx+ 2 dV ’ dx- 2 dV’ 

while the new covariant time derivative operator writes 

d 1 — i d I + i d 
dt 2 2 dt- 


Let us write the stationary action principle in terms of the Lagrange function, as written in Eq. (A13) 


pt2 

SS = 6 L{x,x+,x-,t) dt = 0. 

Ju 


It becomes 


^dL ^ dL dL ^ 

-TT- OX + TT. - 0X+ + TT- - OX- I dt = {). 

ox OXJf- OX— 


Since (5i;+ = d{6x)/dt+ and 5x- = d{Sx)/dt-, it takes the form 


^dL ^ dL 


1 — i d 1-l-z d 
2 2 dt- 


5x\ dt = 0, 


i.e., 




(A15) 

(A16) 

(A17) 

(A18) 

(A19) 

(A20) 

(A21) 


Therefore, after the integration by part, the integral reduces to 



didv) 


Sx dt = 0. 


(A22) 


Finally the Euler-Lagrange equations write 


d dL dL 
dt dV dx' 


(A23) 


They take exactly the form one obtains by writing a stationary action principle based on Eq. ( |A14| ). Moreover, 
once done the transformation d/dt d/dt, this form is nothing but the standard classical form. This result reinforces 
the identification of our tool with a “quantum-covariant” representation, since, as we have shown in previous works 
and as we recall in Sec. this Euler-Lagrange equation can be integrated in the form of a Schrodinger equation. 


APPENDIX B: QUATERNIONIC CALCULUS 

In this appendix, we summarize our notations for the quaternionic algebra and introduce the elementary properties 
of quaternionic arithmetic and analysis which are of use in the present article. 
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1. Definitions and algebraic properties 


A bi-quaternion (p = {(pQ, pi,p 2 ,(p 3 ) is an ordered quadruple of complex numbers. The pi's are the components of 
p. The equality of two quaternions is equivalent to the equality of their corresponding components: p = p H and only 
if pi = with i = 0,1, 2, 3. 


The multiplication of a bi-quaternion by a complex number a, and the addition and substraction of quaternions 
are defined by 

ap = {apo,api, ap 2 , aps ,), 

P + P={Pq+ Pq, P^ + pi,p2 -b p2, p3 + Ps), 

p-p = p+{-l)p. (Bl) 

Addition of quaternions is commutative and associative. The null quaternion p = 0, i.e., the neutral element for 
addition, writes p = (0, 0,0, 0). Multiplication of quaternions by complex numbers is commutative, associative and 
distributive. The quaternionic product pp of two quaternions is itself a quaternion. The product is distributive and 
satisfy, for any complex number a, 

{ap)p = p{ap) = a{pp). (B2) 

Any quaternion can be decomposed as 

p = po + eipi + e 2 p 2 + e3p3, (B3) 

where Hamilton’s imaginary units satisfy the associative but non-commutative algebra 

3 

GiGj — Sij ^ ^ Gij]^€]^, i,j — 1, 2, 3 (^4) 

k=l 

where Cijk is the usual completely antisymmetric three-index tensor with £123 = 1. From these rules, it is easy to 
establish that the product of two arbitrary quaternions is 

pp = {popo - Plpl - P 2 P 2 - p3p3, Popl + PlpO + p2p3 - p3p2, 

P 0 P 2 + P 2 P 0 + p3pl - Plp3, Pop3 + p3p0 + P 1 P 2 - p 2 pl)- (B5) 

This product is not in general commutative, but is associative. 

The adjoint of a quaternion, also named quaternion conjugate or Hamilton conjugate, is defined as 

p^p={p0,-pl,-p2,-p3)- (B6) 


One also defines the scalar product of quaternions 

p-p = ^iPP + PP) = ^iPP + PP) = PoPo + plpl + P 2 P 2 + p3p3- 

The norm of a quaternion is therefore 



(B7) 

(B8) 


The norm and scalar product are numbers, in general complex, the norm being of course the self-scalar product, 
but also the self-quaternionic product. The norm of a product is the product of the norms. When pp = 0, p is said 
to be null or singular. If pp = 1, (/) is unimodular. 


The reciprocal or inverse of any non null quaternion p is defined as 



It possesses the properties: 

pp~^ = p~^p = 1 . 


(B9) 


(BIO) 


The latter property is a key tool of our derivation of the bi-quaternionic Klein-Gordon equation in Sec. fV C 
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2. Quaternions and Conway matrices 


A connection between quaternions and the Dirac-Eddington matrices has first been brought to light by Conway 
1^,^. It gives rise to the definition of a set of sixteen Conway matrices, a subsample of which is used in Sec. ^ to 
derive Dirac’s equation from Klein-Gordon’s. We give below the main steps of their construction. 


Let us consider any two quaternions a and b, and the transformation 

(j) ^ Ip = a4>b, 

which is linear for the elements {(pQ,(pi,(p2i4''i) can therefore be written in matrix form 

tp = M{a, b)(p, 


(Bll) 

(B12) 


where M is a 4x4 matrix and (p and ip are 4x1 column matrices. The notation indicates that M{a,b) is determined 
by a and b. 


If we now consider a couple of linear transformations and their resultant 

%p = a<pb, X = c-ipd, X = cacpbd, 
which can be written in matrix form 

Ip = M{a, b)(p, y = M(c, d)' 0 , y = M(c, d)M (a, b)(p, 

therefore allowing us to write 

M (c, d)M (a, b) = M (ca, bd) 

In particular, for a = c and b = d, Eq. ( |B15 ) becomes 

M{a,bf = M{a'^,b‘^). 


(B13) 

(B14) 

(B15) 

(B16) 


If we substitute a = & = — 1 in Eq. (Bll), we obtain t he id entical transformation and are thus allowed to write 
M(—1, —1) = I, which is the unit matrix. Hence, by Eq. ( B16 ), we have the following theorem: “If a and b are any 
two quaternions satisfying 


a? = b^ = — 1 , 

then, M{a,b)‘^ = /. In other words, the matrix generated by a and & is a square root of unity.” 
Conway proposes the suggestive notation 

M(a, b) = a( )b. 


(B17) 


(B18) 


Now, any of the four quaternions 61 , 62 , 63,1 satisfies Eq. ( B17 ). We have therefore the following sixteen matrices, 
each a square root of unity: 

ei( )6i 6i( )62 ei( )63 ei( )i 

e2( )6i 62) )62 62) )e3 e2( )i 

esi )ei 63) )62 63) )63 63( )i 

*( )ei )62 i( )e3 )i 

It is easy to calculate explicitly these matrices. The whole set of the sixteen Conway matrices can be found in, e.g., 
Ref. [^. We only give below the four matrices used in Sec. ^ 


esi )e2 = 


/ 

0 

0 

0 



/ 

0 

0 

0 

* \ 

0 

0 

-1 

0 

ei( )i = 

0 

0 

—i 

0 1 


0 

-1 

0 

0 


0 

i 

0 

0 

V 

-1 

0 

0 

0 


V 

—i 

0 

0 

0/ 


esi )ei = 


/ 

0 

0 

-1 


/ 

1 

0 

0 

0 \ 


0 

0 

0 

1 1 , . 


0 

1 

0 

0 1 


-1 

0 

0 

0 63) )e3 = 


0 

0 

-1 

0 

V 

0 

1 

0 

0/ 

V 

0 

0 

0 

-1/ 
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3. Quaternionic derivation 


The careful reader will have noticed that, in the whole article, we have avoided to define functions of a quaternionic 
variable which we should have been led to derive with respect to this variable. The reason is simply that, owing to 
the non-commutativity of the quaternionic product, the functions of a quaternionic variable which possess derivatives, 
i.e., for which the right and left derivatives are equal, are only the constant and linear functions. In fact, this is a 
general property shared by every number system with a non-commutative product 


However, the bi-quaternionic wave-function ip, which is defined as a held in the four-dimensional space-time, i.e., 
of which the quaternionic components, ('i/^Q, V'l, V' 2 )' 03 ) of each space-time component ip'^, are complex functions of 
the coordinates X^, can be derived with respect to these coordinates. For instance, d^ip'^ writes 






dip'i 


dX>^ 
3^2 


dXi^ 


62 


dX>^ 


{X°,X\X^,X^) + e3^{X^,X\X^X^). 


(B19) 


We are therefore allowed to apply to the components of this function all the ordinary complex derivation rules with 
respect to the coordinates. 


On the contrary, any relation implying an application to a quaternionic function, other than the constant and linear 
functions, of ordinary derivation rules with respect to the quaternionic variable itself (see, e.g., Eq. (2.52), which is 
however perfectly accurate for a complex ip) would be uncorrect. It is easy to verify that our calculations leading to 
the quaternionic Klein-Gordon and Dirac equations are all along in perfect agreement with this constraint. 
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